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Abstract 

Chronic  kidney  disease  causes  a  slow  loss  of  kidney  function  over  time  and  can  even¬ 
tually  lead  to  End  Stage  Renal  Disease,  where  a  patient  must  undergo  dialysis  to  remove 
fluids  and  wastes  from  the  body.  These  patients  also  suffer  from  a  lack  of  the  hormone 
erythropoietin  (EPO),  produced  in  the  kidneys,  that  stimulates  red  blood  cell  (RBC) 
production.  Without  intervention,  patients  suffer  from  anemia.  Patients  are  treated 
with  both  EPO  and  iron  in  order  to  stimulate  RBC  production. 

We  develop  a  partial  differential  equation  model  for  RBC  dynamics  using  two  struc¬ 
ture  variables,  one  for  age  and  one  for  cellular  iron  endowment.  We  couple  this  with 
a  set  of  ordinary  differential  equations  modeling  iron  dynamics.  We  take  into  account 
the  effects  of  both  inflammation  and  neocytolysis,  which  are  known  to  affect  patients 
undergoing  treatment. 
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1  Introduction 


It  is  estimated  that  31  million  Americans  have  chronic  kidney  disease  (CKD).  Among  those, 
approximately  330  thousand  were  classified  as  being  in  End-Stage  Renal  Disease  (ESRD)  and 
required  dialysis  [42]  to  remove  wastes  and  fluids  from  the  blood.  For  the  purposes  of  this 
study,  we  consider  only  hemodialysis,  where  a  patient’s  blood  is  exposed  to  a  semipermeable 
membrane  outside  of  the  body. 

In  addition  to  regulating  blood  pressure  and  filtering  waste  products  from  blood,  kidneys 
produce  a  hormone  called  erythropoietin  (EPO)  that  is  the  major  regulator  of  erythropoiesis, 
or  red  blood  cell  production.  EPO  level  is  normally  controlled  by  a  negative  feedback  mecha¬ 
nism  in  the  kidneys,  but  patients  in  ESRD  do  not  produce  sufficient  levels  of  EPO  to  maintain 
blood  hemoglobin  concentration.  Hemoglobin  is  the  protein  that  gives  red  blood  cells  the  abil¬ 
ity  to  carry  oxygen.  Patients  with  low  hemoglobin  concentration  may  present  symptoms  of 
anemia,  such  as  decreased  cardiac  function,  fatigue,  and  decreased  cognitive  function. 

In  order  to  prevent  anemia,  patients  typically  receive  recombinant  human  EPO 
(rHuEPO)  intravenously  to  stimulate  red  blood  cell  production.  However,  treatment  is  far 
from  perfect.  In  2006,  only  half  of  dialysis  patients  had  a  mean  monthly  hemoglobin  greater 
than  11  grams  per  deciliter  [42],  the  desired  minimum  level  set  by  the  National  Kidney  Foun¬ 
dation  [29]. 

Iron  is  required  to  produce  hemoglobin,  and  iron  deficiency  can  be  an  issue  among  patients 
receiving  rHuEPO  therapy.  Oral  iron  supplementation  is  often  ineffective,  so  intravenous  iron 
supplementation  has  become  a  mainstay  in  many  patients  undergoing  rHuEPO  therapy  [25]. 

Iron  availability  is  negatively  affected  by  inflammation  level  in  the  body.  Most  patients 
with  CKD  have  elevated  levels  of  inflammation  due  to  CKD  and  the  presence  of  other  medical 
issues  (e.g.,  diabetes,  hypertension,  etc.)  [26]. 

We  develop  a  model  for  red  blood  cell  populations  using  two  structure  variables,  one  for 
age-structuring  and  one  to  account  for  a  cell’s  iron  endowment.  We  couple  these  partial 
differential  equations  with  a  system  of  ordinary  differential  equations  that  models  the  iron 
cycle. 


2  Previous  models 

The  process  of  erythropoiesis  has  been  modeled  in  several  physiological  scenarios.  In  [31], 
rHuEPO  therapy  is  considered  in  healthy  individuals.  This  model  incorporates  the  negative 
feedback  to  endogenous  EPO  production.  EPO  is  assumed  to  be  cleared  using  Michaelis- 
Menten  dynamics.  A  similar  model  was  used  to  fit  data  in  rats  [46].  Both  of  these  models  use 
delay  instead  of  age-structured  modeling. 

Both  [6]  and  [7]  use  age-structured  models,  as  does  the  model  described  in  [27],  which 
assumes  that  the  oldest  mature  erythrocytes  will  be  destroyed,  yielding  a  moving  boundary 
condition.  In  [2] ,  EPO  is  assumed  to  accelerate  maturation  of  cells  undergoing  erythropoiesis. 
Additionally,  EPO  is  assumed  to  be  consumed  during  the  process  of  erythropoiesis. 

The  model  presented  here  is  a  significant  departure  from  these  models  in  that  it  incorpo¬ 
rates  iron.  Additionally,  we  have  accounted  for  the  effects  of  inflammation  and  neocytolysis, 
a  down-regulating  mechanism  by  which  red  blood  cells  are  selectively  destroyed  when  blood 
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oxygen  level  rises. 

The  authors  have  previously  presented  a  more  simplified  version  of  this  model  in  [5,  8], 
including  simulations.  The  simplified  model  is  significantly  less  sophisticated;  assumptions 
were  made  there  so  that  simulations  would  be  more  tractable. 


3  Overview  of  Red  Blood  Cell  Production 

Erythrocytes  are  produced  primarily  from  stem  cells  in  bone  marrow.  In  the  presence  of 
certain  hormones,  stem  cells  divide  asymmetrically,  producing  a  committed  colony-forming- 
unit  (CFU)  while  maintaining  the  population  of  stem  cells.  Erythrocyte  lineage  continues 
as  depicted  in  Figure  1:  erythroid  burst-forming  unit  (BFU-E),  erythroid  colony- forming- 
unit  (CFU-E),  proerythrocyte,  basophilic  erythrocyte,  polychromatic  erythroblast,  orthochro- 
matic  erythroblast,  reticulocyte,  and  erythrocyte. 
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Figure  I:  Erythropoiesis  cell  lineage. 


Hemoglobin  is  synthesized  beginning  in  the  basophilic  erythrocyte  stage,  with  the  majority 
of  synthesis  occurring  in  the  polychromatic  erythroblast  stage.  When  the  nucleus  is  extruded 
from  the  cell,  the  cell  is  named  a  reticulocyte.  Little  hemoglobin  synthesis  happens  at  the 
reticulocyte  stage,  and  synthesis  is  completely  absent  in  mature  erythrocytes  [20].  Reticulo¬ 
cytes  begin  to  lose  the  adhesive  proteins  that  hold  them  in  the  bone  marrow.  They  decrease 
in  size  and  begin  to  circulate  in  the  blood.  In  healthy  individuals,  erythrocyte  life  span  is 
approximately  120  days,  at  which  time  aging  erythrocytes  are  enveloped  by  macrophages  in 
the  spleen. 
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In  a  healthy  individual,  an  increase  in  blood  oxygen  level  is  detected  by  the  kidneys,  which 
in  turn  decrease  production  of  EPO.  Down- regulation  of  the  RBC  production  is  achieved  by 
the  increased  death  rate  of  RBC  progenitors  due  to  this  decrease  in  EPO.  However,  it  takes 
18  to  21  days  for  cells  to  mature  and  be  released  into  circulation,  so  this  down-regulation  has 
a  long  delay  before  it  actually  affects  blood  oxygen  level.  Additionally,  down-regulation  would 
only  cause  the  circulating  red  cell  mass  to  drop  by  about  1%  per  day  [32], 

Neocytolysis  has  been  observed  in  astronauts  entering  space  [4,  43]  and  in  high  altitude 
dwellers  who  descend  to  sea  level  [28,  34];  all  have  higher  aggregate  RBC  mass  than  needed. 
Their  RBC  mass  was  decreased  by  10-15%  in  a  few  days,  which  could  not  be  explained  by 
simple  down-regulation  of  RBC  production. 

Neocytolysis,  believed  to  be  caused  by  a  drop  in  EPO  level,  is  a  physiological  process 
that  aids  in  this  control  of  the  red  cell  mass  by  causing  the  selective  destruction  of  young 
circulating  RBCs  [11,  32],  Since  these  cells  were  circulating  and  contributing  to  blood  oxygen 
level,  their  death  has  a  much  faster  effect  on  blood  oxygen  level  than  the  down-regulation  of 
the  production  cycle.  Neocytolysis  has  been  determined  to  contribute  to  the  anemia  of  renal 
disease  [3,  33],  as  patients  undergoing  therapy  have  constantly  fluctuating  EPO  levels. 

4  RBC  Model 

Cells  in  the  RBC  lineage  divide,  differentiate  and  die  under  the  influence  of  several  hormones. 
Figure  2,  a  modification  of  a  figure  in  [22] ,  depicts  different  aspects  of  the  system  we  will  model. 
This  figure  contains  possible  variables  that  one  might  use  in  formulating  mathematical  models 
for  red  blood  cell  production.  In  particular,  this  figure  suggests  five  cell  classes: 

•  Pi  :  early  BFU-E 

•  P2  :  late  BFU-E,  CFU-E,  proerythroblasts,  and  early  basophilic  erythroblasts 

•  P.i  :  late  basophilic  erythroblasts  and  early  polychromatic  erythroblasts 

•  P4  :  late  polychromatic  erythroblasts,  orthochromatic  erythroblasts,  and  non-circulating 
reticulocytes 

•  P5  :  mature  erythrocytes  and  circulating  reticulocytes. 

We  discuss  each  compartment  individually  in  Sections  4.2  through  4.6,  but  it  should  be 
noted  that  the  original  system  chosen  to  name  the  erythroid  cell  progression  was  based  on 
distinguishing  characteristics  when  samples  were  taken,  stained,  and  viewed  through  a  mi¬ 
croscope.  It  is  not  surprising  that  we  choose  our  five  classes  of  cells  differently  as  we  focus 
on  the  interactions  with  EPO  and  iron  and  the  location  of  the  cells  (in  the  bone  marrow  or 
circulating) . 

The  current  RBC  lineage  model  is  represented  in  Figure  3.  We  introduce  the  model  for 
the  iron  compartments  in  Section  5. 
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Figure  2:  Stages  of  red  blood  cell  production.  We  use  the  following  short-hand  notations: 
BFU-E  =  burst-forming  unit-erythroid,  CFU-E  =  colony-forming  unit-erythroid,  ProEB  = 
proerythroblast,  BasoEB  =  basophilic  erythroblast,  PolyEB  =  polychromatic  erythroblast, 
OrthoEB  =  orthochromatic  erythroblast,  Retie  =  reticulocyte,  RBC  =  red  blood  cell. 
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Figure  3:  Model  schematic.  Note  that  the  states  variables  are  -Pi(t,pi),  P-zit,  /J2),  Pzit-iH 3,7), 
/U4, 7),  Prjt,  P-5, 7)  and  EPO(t).  The  iron  model  is  described  in  Section  5. 
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4.1  Derivation  of  Model  Equations 

We  use  a  model  with  two  structure  variables,  one  for  maturity  level  and  one  for  cellular  iron 
level.  The  structure  variables  //.,  and  7  are  unitless  and  represent  the  maturity  level  and  iron 
state,  respectively,  of  cells.  Thus,  Ps(t,  ^3,7)  is  the  number  of  cells  (in  billions)  in  population 
P3  at  time  t  that  have  maturity  level  p, 3  and  iron  state  7.  Each  class  is  assumed  to  have  a 
maximum  maturity  level,  (/q)/. 

Throughout,  we  use  ft  to  represent  a  birth/proliferation  rate,  8  to  represent  a  death  rate, 
and  p  to  represent  other  rates  (maturation  rates,  rates  of  treatment,  etc.). 

We  will  assume  that  the  iron  state  structure  variable  7  varies  from  0  to  7/  =  2,  where 
7 fuii  —  -j-  =  1  represents  the  amount  of  iron  the  “typical”  RBC  contains  in  a  healthy  individual 
(this  quantity  is  to  be  estimated  later).  We  will  assume  that  when  a  cell  divides,  each  daughter 
cell  inherits  half  of  the  parent  cell’s  iron.  Therefore,  when  a  cell  with  iron  state  7  divides,  it 
leaves  iron  state  7  and  two  cells  enter  iron  state  | .  Similarly,  for  each  cell  in  iron  class  2 7  that 
divides,  two  cells  enter  iron  class  7. 

We  derive  equations  for  those  classes  that  incorporate  iron  by  considering,  at  some  time 
t ,  an  incremental  “area”  of  size  Ap  by  A7,  with  0  <  7  <  =  1,  as  in  Figure  4. 


iron  state 

♦ 


hemoglobinization 
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y  +  Ay 
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maturation 


_ I  | _ ^  maturity 

(.1  (j.  +  Ajli  level 

Figure  4:  Maturation  and  hemoglobinization  processes. 


We  treat  this  region  as  a  arbitrary  elemental  compartment  and  consider  the  time  rate  of 
change  of  the  cell  population  in  the  compartment. 
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q  n+^i  rn+A^ 

di  L  ./, 


P(t,  £,  Qd^d<^  =  2  (rate  entering  iron  class  7  from  class  27) 
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where  0</i</i/,/3is  the  proliferation  rate,  5  is  the  death  rate,  p  is  the  maturation  rate, 
and  h  is  the  hemoglobinization  rate  (equivalent  to  the  rate  of  iron  uptake),  which  generally 
depends  on  the  plasma  iron  level  Fepi. 

We  rearrange  these  to  obtain 


o  />7+A7  pn+An  /»7+A7  p/j,+A/i 

/  P(f ,  7  C)«  =  2  /  /  /5P(t,  7  2C)« 

J  /i  ^7  J 11 

/*7+A7  rn+An 


dt 


n/l+lXll 

/  |d4  f>|P(/.7C)^/C 

'7  Jfi 

/7+A7 

P  [P(t,  p  +  Aju,  C)  -  P(i,  p,  C)] 


'7 


p/x+A/z 


h  [ P(t ,  7  +  A7)  -  P(i,  £,  7)]  df , 
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When  we  allow  A/i  — >•  0  and  A7  — >•  0,  we  obtain 

— P(t,  /i,  7)  =  2f3P(t,  /i,  27)  -  [/3  +  <5]P(t,  /i,  7)  -  Q^Pp{t,  /7  7) 

o 

-  QZhP(^  P  7),  0  <  7  <  ^  =  1,  (1) 

the  general  form  of  the  equations  governing  cells  in  classes  P3  through  P5.  Note  that  the  rates 
/3,  5,  p  and  h  could  be  functions  that  depend  on  p,  7  or  other  states  in  the  system. 

By  our  assumptions,  the  maximum  iron  state  is  7 7  =  2,  so  cell  division  cannot  yield  a  cell 
with  iron  state  7  >  A  =  1.  Thus,  for  cells  with  iron  states  greater  than  A  =  1,  we  omit  the 
birth  rate  term  associated  with  iron  state  27,  as  below: 

r\  r\  O 

Qj.p&  /b  7)  =  -\P  +  S\P(t,  p,  7)  -  Q^pPi1’  Pi  7)  -  Q^hP^i  p  7),  y  =  1  <  7  <  7/  =  2. 

o/i  '  '  “  (2) 


Thus  our  generic  balance  laws  for  the  P3  through  P5  compartments  can  be  written  (using 
the  characteristic  function  y  for  the  interval  [0,1],  i.e. ,  X[o,i]  (7)  —  1  when  7  e  [0,1]  and 
y(q)  =  0  otherwise)  as 

d  d 

— P(t,/i,7)  =  X[0;i]  (/y)2f3P{t,  p,  27)  -  [P-f  h]P(t,/i,7)  -  —pP{t,p,  7) 

-  ^h-Pfa  Pi  7)>  0  <  7  <  2.  (3) 

Cells  in  classes  Pi(t, /ii)  and  P>(b pp  do  not  incorporate  iron,  so  they  are  functions  of 
only  time  and  maturity  level.  The  general  form  of  the  partial  differential  equations  governing 
these  classes,  derived  in  [5],  is 

r\  O 

— P(t,  /i)  =  [/?  -  5\P(t,  p)  -  Qj^pPfa  p)i  (4) 

where  (3  is  the  proliferation  rate,  S  is  the  death  rate,  and  p  is  the  maturation  rate. 


4.2  Class 

Class  Pi  consists  of  early  BFU-E.  These  are  the  most  immature  cells  that  are  committed  to 

the  erythroid  lineage.  We  make  the  following  assumptions  about  cells  in  this  class. 

1.  The  maturation  rate  is  constant,  p\. 

2.  BFU-E  differentiate  into  CFU-E  in  approximately  seven  days  [35,  38].  BFU-E  begin 
expressing  EPO  receptors  (EPORs),  which  means  that,  in  time,  they  do  become  influ¬ 
enced  by  EPO.  Hence,  we  allow  for  cells  to  reside  in  class  Pi  for  three  days  (during 
which  they  are  not  influenced  by  EPO),  while  late  BFU-E  are  in  class  P2  (where  they 
are  under  the  influence  of  EPO).  Therefore  we  choose  the  maximum  maturity  level  to 
be  (pi )f  =  3. 

3.  The  majority  of  cell  proliferation  happens  in  later  classes  [22,  23],  so  we  set  birth  rate 
equal  to  death  rate  for  cells  in  this  class. 

Hence,  for  class  Pi,  equation  (4)  becomes 

d  d 

— Pi(t,/ii)  = -pi— Pi(t,/ii)  (5) 

with  initial  condition 

Pi(o,pi)  =  prVi) 

and  boundary  condition  at  pi  =  0 

P1(t,0)  =  Ph1dy(t). 


4.3  Class  P2(C/o) 


Class  P2  consists  of  late  BFU-E,  CFU-E,  proerythroblasts,  and  early  basophilic  erythroblasts. 
We  make  the  following  assumptions  about  cells  in  class  P2. 


1.  The  maturation  rate  is  constant,  p2. 

2.  Cells  reside  in  this  class  for  twelve  days  [17,  35],  and  therefore  (p2 )/  =  12. 

3.  During  the  time  cells  reside  in  class  P2,  they  undergo  approximately  6  cell  divisions 
[17,  18,  19,  22],  We  assume  that  the  birth  rate  is  a  constant,  /32. 


4.  Cells  in  this  class  express  EPORs  [23,  37],  and  this  interaction  of  EPO  with  EPOR  is  the 
most  important  control  point  for  erythropoiesis  [17,  36,  47].  Cells  at  this  stage  depend 
absolutely  on  EPO  for  their  survival  and  will  undergo  apoptosis  in  its  absence  [23,  38]. 
Even  in  healthy  individuals,  complete  survival  of  progenitors  would  require  an  EPO  level 
much  higher  than  normal  level;  the  normal  production  rate  of  RBCs  represents  survival 
of  only  a  minority  of  progenitor  cells  [22],  Therefore  we  choose  the  death  rate  to  be  a 
decreasing  function  of  the  EPO  level: 


62(EPO) 


zmin 

°2 


k§9 


k§2  /  \  ks<2 

+ (ppoJ 


cg2 


+  (55 
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Thus,  using  equation  (4),  the  state  equation  for  class  P2  is  given  by 

r\  r\ 

— p2(t,  1*2)  =  [P 2  -  52(EPO)]P2(t ,  /i2)  -  p2—P2(t, /i2),  (6) 

with  initial  condition 

p2(o,/i2)  =  prt(/i2) 

and  boundary  condition  at  {p\)f  =  3,  /i2  =  0 

P2(t,0)  =  Pi(*,3). 


4.4  Class  P3(CP3,7) 

Class  P3  consists  of  late  basophilic  erythroblasts  and  early  polychromatic  erythroblasts.  We 
assume  the  following. 


1.  The  maturation  rate  is  constant,  p3. 

2.  Cells  reside  in  this  class  for  1  day  [17,  35],  and  therefore  (p3)f  =  1. 

3.  Cells  in  this  stage  continue  to  proliferate  [17],  with  a  constant  rate  of  fi3  =  In  2  days-1, 
which  corresponds  to  one  cell  division  per  day. 


4.  Cells  in  class  P3  continue  to  express  EPORs,  but  the  level  of  expression  declines  sig¬ 
nificantly  as  cells  mature  through  this  class  [12].  Hence,  cells  in  this  class  become  less 
dependent  on  EPO  for  survival.  We  model  this  as 


53(p3,EPO)=Srx(P3) 


CSS 


kso 


C53 


EPO 


where  5™ax ( p3 )  is  a  decreasing  function,  which  we  assume  to  be  affine  for  this  model 
given  by 


Srx(P3)  =  -m3p3  +  b3,  where  m3  >  0,  b3-m3-(p 3)/ =  0. 


The  function  h3(/i3,  EPO),  for  a  fixed  maturity  level  //3,  is  a  decreasing  sigmoid  function 
of  EPO.  That  is,  at  a  fixed  maturity  level,  as  EPO  increases,  the  death  rate  decreases.  If 
we  consider  instead  a  fixed  EPO  level,  then  increasing  maturation  level  causes  a  decrease 
in  the  maximum  death  rate.  That  is,  low  EPO  affects  younger  cells  in  this  class  more 
than  it  affects  more  mature  cells.  Both  the  fixed  maturity  level  and  fixed  EPO  level 
phenomena  are  depicted  in  Figure  5. 

5.  Cells  in  class  P3  express  transferrin  receptors  (Tfr)  and  begin  the  process  of  taking  in  iron 
and  synthesizing  hemoglobin  [17,  22],  The  rate  of  hemoglobinization  (or  equivalently 
iron  uptake)  is  a  function  of  iron  level  Fepi  (in  the  blood  plasma)  and  iron  state  7  (of  a 
cell).  A  full  description  of  the  hemoglobinization  rate,  h3(Fepi,  7),  along  with  definition 
of  the  parameter  Fepi  appears  in  Section  4.7. 
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Class  P3  death  rate  versus  maturity  level 


Class  P3  death  rate  versus  EPO  level 


maturity  level,  jli3 


Figure  5:  Death  rate  in  class  P3  is  a  function  of  both  maturity  level,  p 3,  and  EPO  level.  For 
a  fixed  EPO  level,  death  rate  decreases  as  maturity  level  increases.  For  a  fixed  maturity  level, 
death  rate  is  larger  when  EPO  is  small.  These  plots  were  generated  with  parameter  values 
m3  =  0.2,  63  =  0.2,  cg3  =  3000  and  ks3  =  20. 


Hence,  for  class  P3,  equation  (3)  becomes 


d  d 

— P3(t,/i3,7)  =  X[o,i]  (7)2/33-P3(£,/hn27)  -  [f33  +  S3(p3,  EPO)]P3(t,  p3,^)  -  p3—P3(t,  p3, 7) 


d_ 

dj 


h3(Fepi,^)P3(t,  /x3, 7) 


0  <  7  <  2, 


(7) 


with  initial  condition 

and  boundary  conditions  at  (/x2) /  =  12,  p3  =  0 


and 


Pzit,  0,7) 


P2(f,12),  7  =  0 
0,  7^0 


P3(t,  /i3, 0) 


P^c  12),  p3  —  0 

P3dy,1(t,p3),  p3^0. 


4.5  Class  Pi(t,  ^4,7) 

Class  P4  contains  late  polychromatic  erythroblasts,  orthochromatic  erythroblasts,  and  non¬ 
circulating  reticulocytes,  the  last  stages  of  erythroid  cells  residing  in  the  bone  marrow.  We 
make  the  following  assumptions  about  cells  in  class  P4. 
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1.  The  maturation  rate  is  constant,  p4 . 

2.  Cells  reside  in  this  class  for  2  days  [17,  35],  and  therefore  (p4)/  =  2. 

3.  Cells  in  this  class  have  stopped  proliferating  [17,  22],  Thus,  /34  =  0. 

4.  It  is  well-documented  that  red  blood  cells  have  a  shorter  life  span  in  individuals  with 
iron  deficiency  [10,  13,  14,  44],  Hence,  we  assume  that  any  RBCs  in  class  P4  that  are 
severely  iron  deficient  have  increased  mortality. 

We  begin  by  defining  c4|  (p4)  as  an  increasing  function  of  /i4.  For  simplicity,  we  choose 
an  affine  function, 


cs4(m)  =  m4/i4  +  64,  m4  >0,  64  >  0, 


with  the  added  restriction  that  0  <  Q,  (/i4)  <  (/r4)/  for  all  /i4. 


We  assume  now  that  death  rate  depends  on  maturity  level  and  cellular  iron  state,  as 
below  : 


<540*4,7 )  =  (<ST“ 


+  <5f“ 


Notice  that  when  7  is  relatively  small,  54(/i4, 7)  is  close  to  5™ax,  and  when  7  is  relatively 
large,  S4 ( /i4 , 7)  is  close  to  However,  the  “relativity”  is  affected  by  c<54(/U4)  when 

7  <  7crit,4-  As  cells  mature  through  class  P4  (i.e.,  as  p4  increases),  then  C74(p4)  increases, 
which  means  that  the  death  rate  for  a  given  iron  level  7  also  increases  (we  choose 
k§4  >  1).  Hence,  a  cell  having  only  a  little  iron  when  it  arrives  in  class  P4  is  less  likely 
to  die  than  a  cell  with  that  same  level  of  iron  that  is  about  to  mature  out  of  class  P4. 
Figure  6  depicts  fi4(/r4, 7)  over  varying  p4  and  7. 


5.  Cells  in  class  P4  continue  to  collect  iron  and  synthesize  hemoglobin  [17,  22]  at  a  rate 
h4(Pep/,7),  explained  further  in  Section  4.7. 


Thus,  using  equation  (3),  the  state  equation  is  given  by 

d  d 

— P4(t,  p4, 7)  =  -  <M/W  h4, 7)  -  /N,  7) 


d_ 

dj 


hi(Fevi,  7)P4(t,  /i4, 7) 


(8) 


with  initial  condition 


P4(0,^4,7)  =  Prt(^4,7) 
and  boundary  conditions  at  (^3)/  =  1,  /r4  =  0 

P*(t,  0,7)  =  P3(t,l,7) 


and 


P4(t,/i4, 0) 


P3(t,  1,0),  yU4  =  0 
P4My’7(t,^4),  jU4  7^  0. 
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Class  P4  death  rate  versus  maturity  level  Class  P4  death  rate  versus  iron  state 


Figure  6:  Death  rate  in  class  P\  is  a  function  of  both  maturity  level,  /.14.  and  cellular  iron 
state,  7.  For  a  fixed  maturity  level  p 4,  death  rate  decreases  as  iron  state  7  increases.  That 
is,  for  a  given  maturity  level,  cells  with  little  iron  are  more  likely  to  die.  For  a  iron  state 
7  <  7cr*/,,4 ,  death  rate  is  increases  as  maturity  level  increases.  These  plots  were  generated  with 
parameter  values  S'™1  =  0.1,  5™lx  =  0.5,  7714  =  0.38,  64  =  0,  and  Ay,  =  12. 


4.6  Class  P5(A,/i5,7) 

Cells  in  class  P5  are  mature  erythrocytes  and  circulating  reticulocytes.  For  cells  in  this  class, 
we  make  the  following  assumptions. 

1.  The  maturation  rate  is  constant,  p5. 

2.  The  average  life  span  of  a  red  blood  cell  in  a  patient  with  chronic  kidney  disease  is  about 
70  days,  significantly  shorter  than  for  healthy  persons  [24,  39,  41].  We  set  the  maximum 
maturity  level  to  be  (^5)/  =  100  days-1,  but  the  death  rate  is  set  such  that  virtually  no 
cells  reach  the  maximum  maturity  level. 

3.  As  they  have  no  nuclei,  cells  in  this  stage  do  not  proliferate  [17],  and  therefore  (35  =  0. 

4.  Death  rate  in  this  class  is  assumed  to  have  four  components:  death  due  to  low  cellular 
iron,  death  due  to  aging,  death  due  to  neocytolysis,  and  death  due  to  blood  loss. 


(a)  Death  due  to  low  cellular  iron.  As  in  the  previous  class,  cells  that  lack  a  full 
complement  of  iron  have  increased  mortality.  As  in  class  P4,  we  have 

c<55.7(h5)  =  m5p5  +  65,  m5>  0, 
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and 


^5,7  (/^5)  T) 


+  5 


min 
5,7  ‘ 


(b)  Death  due  to  aging.  Senescent  (aged)  erythrocytes  are  enveloped  and  destroyed  in 
the  spleen  [17],  so  death  rate  in  this  class  is  also  a  function  of  maturity  level.  This 
process  occurs  independent  of  cellular  iron  state: 


zmax  rmm 

5  5 


I  rmm 


(c)  Death  due  to  neocytolysis.  Neocytolysis  targets  RBCs  between  14-21  days  old  (i.e. , 
for  maturity  level  14  <  /i5  <  21)  [4,  32],  These  cells  have  increased  death  rate  when 
EPO  level  drops,  but  it  is  not  yet  known  whether  the  mechanisms  that  cause  this 
increase  in  death  rate  are  related  to  a  large  drop  in  EPO  level  (i.e.,  a  large  negative 
rate  of  change  of  EPO  level)  and/or  simply  a  low  EPO  level.  We  account  for  both 
in  this  model,  with  tuning  parameters  ay  and  a 2  to  account  for  the  relative  effects. 


We  begin  with  the  assumption  that  neocytolysis  is  a  response  to  low  EPO  level. 
We  choose  a  decreasing  sigmoid  function  for  the  maximum  death  rate: 

uneo 

( rneo\ 

_ ) _ 

.  j uneo  .  .  j uneo  • 

c^°) **  +  (EPO)  *5 

Then  the  death  rate  is  given  by 

£5,ne0>1(/W  EPO)  =  5™*oA(EPO) 

the  sum  of  smoothed  heaviside  functions.  Note  that  when  kh,i  is  chosen  large 
enough,  the  death  rate  is  zero  when  <  14  and  /J5  >21.  Also,  for  maturity  levels 
14  <  ^5  <  21,  when  EPO  level  is  small,  the  maximum  death  rate  is  large.  Figure 
7  depicts  the  output  of  the  death  rate  function  45  ne0)i  (//5,  EPO),  for  both  the  case 
of  a  fixed  maturity  level  with  varying  EPO  and  the  case  of  fixed  EPO  with  varying 
maturity  level. 


1  1 

^  _|_  g— 2fcfi,,i(M5— 15)  1  _|_  g— 24,lfe- 20)  ’ 


zmax 

^5,neo,l 


(EPO)  =  % 


max,! 


For  the  assumption  that  neocytolysis  is  a  response  to  large  negative  rate  of  change 
of  EPO,  we  choose  a  decreasing  function  for  the  maximum  death  rate: 


<5 


max 

5,neo,2 


tepo) 

dt  J 


cmax,  2 
^5,neo 


1 


1 

1  _|_  g.  2A7,.:j(  EPO  pEPO.crit) 


(Note  here  that  we  do  not  use  a  sigmoid  function  (as  we  did  for  8™®*0,i(EPO)) 
because  the  area  of  interest  includes  negative  values  of  ^ EPO .)  The  death  rate, 
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Class  P5  death  rate  due  to  neocytolysis 
(low-EPO  mechanism)  versus  maturity 


maturity  level,  jli 


Class  P5  death  rate  due  to  low-EPO 
neocytolysis  versus  maturity  level 


Figure  7:  Low-EPO  neocytolysis  death  rate  in  class  P5  is  a  function  of  maturity  level  /i5 
and  EPO  level.  In  the  first  plot,  we  fix  EPO  and  show  the  death  rate  over  varying  maturity 
level.  Notice  that  neocytolysis  affects  cells  such  that  14  <  /i5  <  21.  Also,  as  EPO  increases, 
the  effects  of  neocytolysis  diminish.  This  fact  is  demonstrated  again  in  the  second  plot.  To 
produce  these  plots,  we  chose  parameter  values  =  0.2,  c^°  =  2000,  k's^°  =  20  and 


h,  1  =  10. 


given  by 


& 5,neo,2  (  5;  ^  EPO 


=  %S*[±EP0) 


1 


1 


2  _|_  g  2/^2  (^5  15)  2  -|-  g  20) 


is  plotted  in  Figure  8  for  both  the  constant  EPO  rate  with  varying  maturity  case 
and  the  constant  maturity  level  with  varying  EPO  rate  case. 


Then  the  death  rate  due  to  neocytolysis  is  given  by 


^5  ,neo 


h5 j  EPO ,  —EPo\  —  ai55;ne0ji  (/X5,  EPO )  +  0:2^5, neo, 2  ( h 5,  —  EPC)  \  . 


dt 


d 


dt 


(d)  Death  due  to  blood  loss.  We  make  a  separate  death  rate  to  account  for  losses  such 
as  blood  lost  during  blood  draws,  denoted  5c,:ioss(t,  ^5).  The  iron  from  these  losses 
is  not  recycled. 
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Class  P5  death  rate  due  to  neocytolvsis 
(EPO-rate  mechanism)  versus  maturity 


maturity  level,  p 


Class  P5  death  rate  due  to  neocytolvsis 


EPO  rate,  units  of  EPO  per  day 


U 

a 

£2 


Figure  8:  Large  negative  EPO  rate  neocytolysis  death  rate  in  class  P5  is  a  function  of  maturity 
level  /.(r,  and  EPO  rate.  Note  that  neocytolysis  targets  cells  with  maturity  level  14  <  <  21. 

As  EPO  rate  attains  large  negative  values,  the  death  rate  due  to  neocytolysis  increases.  To 
produce  these  plots,  we  chose  parameter  values  <5™nef  =  0-2,  0i,,3  =  0.0075,  pEPO,crit  =  750 
and  kh,2  =  5. 


We  assume  that  these  phenomena  occur  independently  and  that  the  resultant  death  rate 
for  class  P5  is  the  sum  of  these  four  components, 


^5  7 j  EPO ,  —EPO^j  —  ^5)7(^5, 7)  +  (^(/u) 


d 


$5 ,neo  (  h5)  EPO,  ~EPO 


^5 doss 


5.  Only  reticulocytes  at  this  stage  are  still  able  to  collect  iron,  and  at  a  smaller  rate  than 
previous  cell  classes  [9,  15,  17].  As  before,  the  hemoglobinization  rate,  /i5 ( Fepi , 7) ,  for 
cells  in  class  P5  is  discussed  in  Section  4.7. 


Using  equation  (3),  cells  in  class  P5  are  governed  by  the  equation 

^P5(t,  /i5, 7)  =  -S5  ^/i5, 7>  EPO,  ^EPO^j  P5(t,  fi5, 7)  -  P5^ 


d_ 

<9y 


hs(Fevi,"t)P$(t, 


(9) 


with  initial  condition 


Ps{o,Kn)  =  Pta(K,'i) 

and  boundary  conditions  at  (/i4) /  =  2,  /i5  =  0 

^5^,  0,7)  =  Pi(t,  2, 7) 


16 


and 


P5(t,  ^5,0) 


P^t,  2,0),  n 5  =  0 


4.7  Hemoglobinization 

Cells  in  classes  P3,  P4  and  P5  participate  in  iron  uptake  for  the  purposes  of  synthesizing 
hemoglobin.  As  noted  earlier,  we  use  7  as  a  unitless  structure  variable  that  indicates  the  level 
of  iron  in  a  given  cell.  Moreover,  recall  that  0  <  7  <  2,  where  ^fuii  =  1  represents  the  amount 
of  iron  the  “typical”  RBC  contains  in  a  healthy  individual.  Finally  recall  when  a  cell  divides, 
we  assume  that  each  of  the  daughter  cells  receives  half  of  the  parent  cell’s  iron  endowment. 

In  order  to  determine  the  rate  of  iron  moving  from  the  iron  compartment  and  entering 
RBCs,  we  compare  the  rate  of  iron  required  for  cells  to  be  hemoglobinized  at  the  maximum  rate 
with  the  rate  of  iron  available  in  the  blood  plasma  (i.e.,  iron  in  the  blood  plasma  compartment). 
Then  we  can  determine  if  there  is  sufficient  iron  available  to  hemoglobinize  at  the  maximum 
rate  or  if  hemoglobinization  is  occurring  in  an  iron-restricted  fashion. 

We  assume  that  when  there  is  an  abundance  of  iron  available,  cells  can  obtain  a  full  com¬ 
plement  of  iron  in  three  days  [17].  Cells  begin  the  hemoglobinization  process  in  class  P3.  Cells 
reside  in  this  class  for  one  day  and  they  undergo  cell  division  at  a  rate  of  1  division  per  day. 
After  cells  transition  out  of  class  P3,  they  no  longer  undergo  cell  division  (i.e.,  cells  in  classes 
P4  and  P5  do  not  divide).  While  we  model  hemoglobinization  and  proliferation  as  continuous 
processes,  the  following  considerations  in  discrete  time  nodes  indicate  that  the  maximum  rate 
of  iron  uptake  should  be  set  to  fuu  days-1  so  that  a  cell  could  attain  a  full  iron  complement 
in  three  days’  time. 


day  0 

day  1  before  cell  division: 
day  1  after  cell  division 
day  2 
day  3 


7  =  0 

7  =  \lfuii  =  | 

7  =  b fun  =  l 
7  =  \lfull  +  \lfull  =  | 

7  =  §7 fun  +  fuii  —  1 


Cells  in  class  P4  can  uptake  iron  at  this  rate  until  their  iron  state  reaches  7/„«,  at  which 
point  they  stop  taking  in  iron.  Reticulocytes  in  class  P5  (i.e.,  cells  in  class  P5  with  maturity 
level  0  <  /i5  <  2)  are  circulating  in  the  blood.  While  they  are  still  capable  of  iron  uptake 
until  they  reach  iron  state  7 fua  =  1,  the  rate  of  uptake  is  smaller  [17],  so  we  multiply  the  rate 
for  class  P4  by  the  constant  0  <  Ay  <  1.  Therefore  we  define  the  maximum  iron  uptake  rate 
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functions  h™ax{ 7),  h™ax( 7),  and  h™ax(n 5,7),  depicted  in  Figure  9,  by 


&r*(7) 

^(7) 


2 

5 


2  1 

5  1  +  ,max 


f  k5hrx( 7),  0  <  /i5  <  2 

1  0,  /i5  >  2. 


Then  the  rate  of  iron  needed  if  every  cell  becomes  hemoglobinized  at  its  maximum  rate  is 

PFe, needed  =  (max  rate  iron  needed  in  class  P3)  +  (max  rate  iron  needed  in  class  P4) 

+  (max  rate  iron  needed  in  class  P5) 


=  k 


Fe 


Hf  r1  2 

/  -P3(t,/i3,7)d/i3d7 


to  Jo 


+ 


[  [  h?axh)P,(U  IM,  jjdfMdr/ 

'0  Jo 


Hf  r 2 


+ 


Kax{p  5, 7)^5^,  P5,  i)dpdi 


'0  Jo 


where  kpe  is  the  amount  of  iron  in  one  billion  red  blood  cells  with  a  full  iron  endowment. 

We  compare  this  with  the  rate  of  iron  available  for  erythropoiesis  at  time  t.  We  assume 
that  the  iron  available  for  erythropoiesis  is  proportional  to  amount  of  iron  in  the  blood  plasma 
compartment.  Hence,  the  rate  of  iron  availability  is  given  by 


PFe, avail 


^Fe,ef  f  FCpl  (t) , 


Maximum  Iron  Uptake  Rates  for  classes  RBC 


classes  P3,  P 4  and  P5 


Figure  9:  Maximum  iron  uptake  rates  where  0  <  p$  <  2. 
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where  0  <  kFe,eff  <  1  days  1 . 

We  define  the  iron  availability  rate  fraction,  0  <  kpe, avail  <  1,  to  be 

PFe, needed  —  PFe,  avail  i 
pFe, needed  ^  P Fe, avail- 

Then  the  actual  rates  of  iron  uptake  are  given  by 

hi(Fepi, 7)  =  kFe,avaiih™ax( 7),  i  =  3, 4,  5. 


( 

kFe,  avail  =  <  PF^cmari 
\  PFe, needed, 


4.8  EPO  compartment 

While  healthy  individuals  produce  EPO  in  response  to  a  decrease  in  blood  oxygen,  patients 
whose  kidneys  have  only  minimal  function  produce  only  a  small  basal  level  of  EPO  in  the 
kidney  and  liver  [21],  They  receive  rHuEPO  intravenously  during  their  dialysis  treatments. 
In  this  model,  we  assume  the  use  of  a  bioidentical  rHuEPO;  we  will  not  distinguish  between 
rHuEPO  and  endogenous  EPO  with  respect  to  their  action. 

We  assume  the  rate  of  endogenous  EPO  production  in  the  liver  and  kidney  to  be  constant. 
We  will  assume  that  EPO  clearance  is  proportional  to  the  amount  present.  Thus  we  have 

—  EPO  {t)  =  pEPO,endog  +  PEPO,exog(t )  —  - - EPO(t),  (10) 

at  1 1/2 

with  initial  condition 

EPO(  0)  =  £0, 

where  pFpo,endog  is  the  (assumed  constant)  rate  of  endogenous  EPO  production,  pEPO.exog(t) 
is  the  rate  of  exogenous  EPO  provided  during  treatment,  and  ti/2  is  the  half-life  of  EPO. 

5  Iron  Model 

We  now  introduce  the  second  portion  of  our  model,  the  iron  compartments.  Iron  is  required 
to  make  hemoglobin,  the  protein  that  gives  erythrocytes  the  ability  to  carry  oxygen.  It  is 
also  the  protein  that  gives  erythrocytes  their  characteristic  red  color.  If  iron  is  not  available 
during  erythropoiesis,  the  result  is  lighter-colored  (hypochromic)  erythrocytes  with  reduced 
capacity  to  carry  oxygen. 

Control  of  iron  in  the  body  is  a  strictly  regulated  process,  in  part  because  there  is  no 
pathway  for  the  excretion  of  excess  iron  [16,  45].  When  red  blood  cells  age,  they  become 
enveloped  by  macrophages  in  the  spleen.  The  iron  from  their  hemoglobin  is  then  recaptured 
and  sent  to  the  bone  marrow  for  use  in  making  hemoglobin  for  new  erythrocytes.  This 
recycling  process  is  very  efficient  and  is  the  main  source  of  iron  to  erythropoiesis  [40].  In 
much  smaller  quantities,  iron  is  absorbed  from  diet  and  is  transferred  in  and  out  of  storage  in 
the  liver.  The  only  losses  to  the  system  are  from  sweating,  cells  being  shed,  blood  losses,  etc. 

The  major  regulator  of  iron  homeostasis  is  the  hormone  hepcidin,  which  is  produced  in 
the  liver  [16,  30,  45].  Hepcidin’s  role  in  iron  homeostasis  has  only  been  recently  studied;  the 
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(a)  Ferroportin  is  required  for  the  transport 
of  iron  out  of  cells. 


(b)  Hepcidin  is  the  major  regulator  of  iron 
transport  out  of  cells. 


Figure  10:  Iron  regulation  at  a  cellular  level. 


mechanisms  and  pathways  leading  to  the  production  and  action  of  hepcidin  are  only  beginning 
to  be  understood. 

The  protein  ferroportin  is  required  to  transport  iron  out  of  a  cell  and  into  the  plasma. 
Hepcidin  binds  to  ferroportin  and  causes  the  complex  to  be  absorbed  into  the  cell,  effectively 
interrupting  the  transport  of  iron  into  the  blood  plasma,  as  depicted  in  Figure  10.  Hence,  an 
increase  in  hepcidin  level  causes  iron  to  remain  in  cells  and  not  be  released  into  the  blood 
plasma  (where  it  could  be  used  in  RBC  production). 

Hepcidin  production  is  increased  in  response  to  high  iron  level  in  the  blood  plasma  and/or 
liver  and  is  decreased  in  response  to  erythropoietic  activity  [16,  30].  Production  is  also  in¬ 
creased  in  the  presence  of  certain  cytokines  which  are  released  due  to  inflammation  in  the  body 
[16,  45].  It  is  thought  that  this  might  be  a  defense  mechanism  against  foreign  organisms  which 
may  need  iron  to  reproduce.  Since  patients  in  ESRD  commonly  have  other  health  problems, 
they  often  have  higher  than  normal  levels  of  inflammation.  Thus,  they  may  produce  higher 
than  normal  levels  of  hepcidin  and  as  a  result,  even  if  there  is  enough  iron  in  the  body,  it  may 
not  be  available  for  erythropoiesis  because  it  cannot  leave  the  cells  and  enter  the  plasma  [1]. 

We  propose  a  model  of  ordinary  differential  equations  for  the  iron  model  using  the  con¬ 
ceptual  model  given  in  Figure  11. 

5.1  Hepcidin 

We  choose  H(t )  to  represent  plasma  hepcidin  level  at  time  f,  with 

0  <  H(t)  <  1. 

A  hepcidin  level  of  0  indicates  little  or  no  hepcidin  present.  In  this  situation,  iron  can  pass 
unencumbered  across  a  cell  membrane  and  enter  the  iron  in  the  blood  plasma  compartment. 
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\ 


Figure  11:  Iron  Model.  Note  that  the  states  variables  are  H ,  Feu,  Fepi  and  Fere. 


A  hepcidin  level  of  1  indicates  a  maximum  hepcidin  level;  when  hepcidin  is  at  its  maximum 
level,  iron  cannot  cross  a  cell  membrane  to  enter  the  blood  plasma.  In  order  to  ensure  that 
0  <  H  <  1,  we  use  characteristic  functions  that  are  “turned  on”  when  a  condition  is  true  and 
“off’  when  it  is  false.  For  example, 

m  /  r  H$)  <  r 

=  X{r\H{r)<l}{t)  =  <  Q  >  ^ 

We  assume  the  following: 

1.  Hepcidin  is  produced  in  response  to  high  levels  of  iron;  how  this  is  measured  by  the  body 
is  unclear.  We  assume  that  production  is  increased  in  response  to  high  levels  in  both 
the  liver  and  in  the  plasma.  The  relative  size  of  these  two  responses  can  be  adjusted  by 
way  of  the  coefficients  kpe,pi  and  kpe,ii  ■ 

rate  of  increase  in  hepcidin  due  to  high  level  of  iron  in  plasma  =  k pe )P; Fepi (t)x{H<i}> 
rate  of  increase  in  hepcidin  due  to  high  level  of  iron  in  liver  =  kFe,uFeu(t)x{H<i}, 

with  0  <  kf-p  pi ,  kpe,ii  <  1.  The  characteristic  functions  in  these  terms  allow  hepcidin  level 
to  increase  to  H  =  1  in  response  to  high  level  of  iron  in  the  plasma  or  liver,  but  do  not 
allow  hepcidin  level  to  increase  above  H  —  1. 

2.  Hepcidin  production  is  increased  in  the  presence  of  inflammation,  but  the  specific  path¬ 
ways  that  cause  this  response  are  still  under  investigation.  We  assume  that  the  input 
I(t)  is  a  measure  of  the  overall  inflammation  state  of  the  body  (a  broad  simplification 
that  should  be  revisited  in  future  model  iterations)  such  that  0  <  I  <  1.  Therefore, 

rate  of  increase  in  hepcidin  due  to  inflammation  =  kinfiI(t)x{H< i}, 
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where  0  <  kinfi  <  1. 

3.  While  hepcidin  is  actually  “consumed”  as  it  is  used  to  prevent  the  action  of  ferroportin, 
we  initially  assume  a  clearance  proportional  the  amount  present. 

Hence,  we  model  hepcidin  level  by 

dH 

-^r(£)  =  kinfiI(t)x{H<i}(t)  +  kFe,piFepi(t)x{H<i}(t)  +  kmiFeH{t)x{H<i}(f)  ~  kHH(t)x{H>o}(t), 

(H) 

where  0  <  kn  <  1. 


5.2  Iron  in  the  reticuloendothelial  system 

The  reticuloendothelial  system  consists  of  macrophages  in  the  spleen  and  other  cells  that 
collect /envelop  old  and  damaged  cells  and  reclaim  their  iron.  This  includes  any  cell  in  the  red 
blood  cell  lineage  that  dies,  as  well  as  liver  cells  that  become  senescent  and/or  damaged. 

Based  on  our  model  for  RBC  dynamics,  the  rate  of  iron  in  to  the  reticuloendothelial  system 
from  RBC  lineage  is 


vFe 


Hf  r1 


53(/x3,  EPO)^P3{t ,  p3, 7)d/x3d7 


'0  JO 


rif  r  2 


+ 


hi(p4,  EPOFP4(t,  /r4, 


1 0  J  0 


clOO  r 


+ 


'0  Jo 


^5,7 (h-5)  T)  T 


d 


$5 ,neo  (  P 5)  EPO,  —EPO 
'  at  J 


'fP5(t,n5,'y)dn5d'y  , 


where  kpe  is  the  amount  of  iron  in  one  billion  RBCs  with  a  “normal”  iron  endowment.  Notice 
that  for  class  P5  we  do  not  recycle  iron  from  deaths  due  to  blood  losses. 

The  liver  serves  as  a  storage  area  for  iron.  Iron  residing  in  liver  cells  that  die  is  recaptured 
and  ends  up  in  the  reticuloendothelial  system.  Hence,  we  assume  that 


rate  of  iron  in  to  the  reticuloendothelial  system  from  the  liver  =  ku^reFeu(t) , 


where  0  <  ku^re  <  1. 

Reticuloendothelial  cells  require  ferroportin  to  transfer  iron  from  within  the  cell  into  the 
blood  plasma.  Hepcidin  negatively  regulates  this  transport;  when  hepcidin  level  is  high,  the 
amount  of  iron  that  can  leave  liver  cells  and  enter  the  blood  plasma  is  low.  We  assume  that 


rate  of  iron  out  of  the  reticuloendothelial  system  in  to  the  blood  plasma  =  kre^pi{l—H{t))Fere{t), 
where  0  <  kre^pi  <  1. 
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Thus,  we  have 


dFere 

dt 


vFe 


Hf  /'! 


Jo  Jo 

•7/  r2 


+ 


S3(/J3,  EPO)^Pz{t,  fiS,  7)^3^7 
^i(/-W  EPO)~fP4(t,  im,  7)0^7 


'0  ^0 
/•7/  />100  r 


+ 


*^5,7 ( h'5 ;  7)  T  ^U,/x C /^j7j )  T  ^5 ,neo  ^/h>>  EPO ,  ^EPO^ 


JO  JO 

T  kn^reE&ii  ^Ve— *-/>/(  f  H)Fere. 


7P5(£,/i5,7)d/i5d7 

(12) 


5.3  Iron  in  the  Liver 

Iron  enters  the  liver  (for  storage)  from  the  blood  plasma  when  the  level  of  iron  in  the  blood 
plasma  is  high  and  leaves  storage  when  the  level  of  iron  in  the  blood  plasma  is  low.  We  model 
these  rates  with  sigmoid  functions  as  follows: 


rate  of  iron  into  liver  from  the  blood  plasma  compartment  =  Kpi^u  ( Fepi(t ))  Fepi(t), 


where 


and 


Epi—tii  ( Fe.pi )  kri 


(Cpl—tli) 


kpl—tli 


(cpi^ii)  pl^li  +  (Fep/) 


jpl—±li 


rate  of  iron  out  of  the  liver  into  the  blood  plasma 


=  Ku^pi  ( Fepi(t ))  (1  -  H(t))Feu(t) 


where 


Ku—tpi  ( Fepi ) 


(cH^pl)ku^1 


(cu^pi)kli^pl  +  (Fepi) 


kli—tpl 


Notice  that  1  —  H  appears  in  the  expression  for  rate  of  iron  into  the  blood  plasma  compartment 
from  the  liver  compartment  because  hepcidin  inhibits  the  action  of  ferroportin  in  allowing  the 
transport  of  iron  out  of  the  liver.  The  functions  kpi^u  ( Fepi )  and  ku^pi  ( Fepi )  are  depicted  in 
Figure  12. 

When  cells  in  the  liver  die,  the  iron  in  those  cells  is  collected  by  cells  in  the  reticuloen¬ 
dothelial  system.  We  model  this  as 


rate  of  iron  out  of  the  liver  in  to  the  reticuloendothelial  system  =  ku^.reFeu(t) . 


Hence,  the  time  rate  of  change  of  iron  in  the  liver  is  given  by 


dFeu 

dt 


Kpi^u  ( Fepi(t ))  Fepi(t)  -  Ku^pi  ( Fepi(t ))  (1  -  H(t))Feu(t)  -  ku^reFeu(t) . 


(13) 
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—  .  —  Kpi^n(Fepi(t)) 
Kn^pi(Fepi(t)) 


^  mg  iron 


Figure  12:  Rate  coefficients  for  iron  in  and  out  of  storage  (i.e.,  iron  into  the  liver  from  the 
blood  plasma  and  iron  into  the  blood  plasma  from  the  liver). 


5.4  Iron  in  the  Blood  Plasma 

Iron  in  the  blood  plasma  is  the  main  source  of  iron  for  erythropoiesis.  Iron  added  to  blood 
plasma  during  intravenous  treatment  is  denoted  Feex(t).  Iron  also  enters  the  blood  plasma 
from  the  reticuloendothelial  system,  the  liver,  and  the  intestines  from  diet.  Since  each  of  these 
requires  ferroportin  to  transfer  the  iron  from  within  the  cell  to  the  plasma,  the  associated  rates 
of  transfer  are  all  subject  to  the  negative  regulatory  effect  of  hepcidin. 

Iron  enters  the  intestines  from  a  patient’s  dietary  intake.  We  define 

Piron,diet  =  maximum  rate  of  iron  coming  from  diet,  mg/day 

Iron  in  intestinal  cells  must  be  transported  into  the  blood  plasma  via  ferroportin.  Hence,  the 
rate  of  iron  into  the  blood  plasma  is  limited  by  the  hepcidin  level.  Thus,  we  assume 

rate  of  iron  into  blood  plasma  from  diet  =  kdiet(  1  —  H) Piron,diet, 

where  0  <  k(jtet  <  1. 

We  define 


Feex(t )  =  rate  of  iron  into  plasma  from  intravenous  treatment 


and 

pFe,pi,ioss (t)  =  rate  of  iron  lost  via  blood  plasma  losses. 

The  rates  of  iron  (i)  into  blood  plasma  from  the  reticuloendothelial  system,  (ii)  into  the 
blood  plasma  from  the  liver,  (iii)  into  the  liver  from  the  blood  plasma,  and  (iv)  into  the  red 
blood  cell  lineage  were  all  described  previously.  Thus,  we  have  the  following  equation  for  iron 
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in  the  blood  plasma: 


dFcpi 

dt 


=  rate  of  iron  into  plasma  from  diet 


+  rate  of  iron  into  plasma  from  reticuloendothelial  system 
+  rate  of  iron  into  plasma  from  liver 

+  rate  of  iron  into  plasma  from  intravenous  treatment 
—  rate  of  iron  into  liver  from  plasma 

—  rate  of  iron  into  RBC  lineage  from  plasma 
—  rate  of  iron  lost  via  blood  plasma  losses 


kdieti^  1 


iron, diet 


+  kre^pi  (1  -  H)Fere(t ) 


+  Ku^pi(  1  —  H)Fen(t ) 

+  Feex(t) 


Fpi—riiiFcpi  ]  Ff'pi  (tj 

k-Fe  [hsiFcpi ,  -|-  h^Fcpi^^  T  h§(Fcpi,  y)] 

PFe,pl,loss\t)  • 


(14) 


6  Summary 

We  have  developed  a  model  for  both  red  blood  cell  and  iron  dynamics  in  patients  in  ESRD 
undergoing  hemodialysis.  In  particular,  our  model  consists  of  the  RBC  precursor  and  cir¬ 
culating  fragmentation  equations  (5)- (9)  with  appropriate  initial  and  boundary  conditions, 
coupled  with  (10)- (14).  We  use  structure  variables  to  account  for  the  age  and  iron  endowment 
of  cells  in  the  red  blood  cell  lineage  in  the  fragmentation  equations  for  P\  through  P5.  The 
model  also  permits  investigation  of  the  effects  of  neocytolysis  and  inflammation  on  red  blood 
cell  dynamics.  The  coupled  fragmentation  and  EPO,  iron,  hepcidin  with  inflammation  in¬ 
put  are  nontrivial  to  solve.  We  are  currently  carrying  out  simulations  on  both  the  continuous 
model  presented  here,  and  a  corresponding  discrete  time  version,  in  verification  and  validation 
studies. 
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